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Abstract 

Background: Vesiviruses in tine family Caliciviridae infect a broad range of animal hosts including mammals, birds, 
fish, amphibians and reptiles. The vesivirus Crol strains were isolated from diseased snakes in the San Diego zoo in 
1978 and reported as the first caliciviruses found in reptiles. The goal of this study was to characterize the Crol 
strain 7800321 that was isolated in cell culture from a rock rattlesnake {Crotalus lepidus) in the original outbreak 

Results: We re-amplified the original virus stock in Vero cells, and determined its full-length genome sequence. The 
Crol genome is 8296 nucleotides (nt) in length and has a typical vesivirus organization, with three open reading 
frames (ORF), ORFl (5643 nt), 0RF2 (2121 nt), and 0RF3 (348 nt) encoding a nonstructural polyprotein, the major 
capsid protein precursor, and a minor structural protein, respectively. Phylogenetic analysis of the full-length 
genome sequence revealed that the Crol virus clustered most closely with the VESV species of the genus Vesivirus, 
but was genetically distinct (82-83% identities with closest strains). 

Conclusions: This is the first description of a full-length genome sequence from a reptile calicivirus (Crol). The 
availability of the Crol genome sequence should facilitate investigation of the molecular mechanisms involved in 
Crol virus evolution and host range. 
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Background 

The family Caliciviridae is a large group of small, non- 
enveloped RNA viruses that includes important human 
and animal pathogens [1]. The family is comprised of 
five genera: Lagovirus, Vesivirus, Nebovirus, Sapovirus 
and Norovirus [2] and two new genera have been pro- 
posed [3]. Despite marked genetic and antigenic diversity, 
caliciviruses share several common features. All have icosa- 
hedral virions with a protein shell containing 180 copies of 
a major capsid protein, VPl [4]. The virions carry a 
positive-sense single-stranded RNA genome of approxi- 
mately 6.4 to 8.5 kb in length. The 3'-end of the calicivirus 
RNA is polyadenylated, and the 5'-end is covalentiy linl<ed 
to a small protein encoded by the virus genome, VPg [5]. 
Calicivirus RNA genomes share similar organization; they 
are comprised either of two or three ORFs. The large 
ORFl encodes the virus nonstructural proteins and is 
expressed from the genomic RNA template. In the gen- 
omes of sapo-, lago-, and neboviruses, the nonstructural 
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region is fused to the gene encoding virus capsid protein, 
VPl. In the genomes of vesi- and noroviruses, the capsid 
protein is encoded by a separate ORF2, located towards the 
3'-end of the virus genome. For all caliciviruses, the capsid 
proteins are produced from an abundant subgenomic RNA 
synthesized during virus replication. The same RNA serves 
as bicistronic template for the expression of a minor capsid 
protein, VP2. The ORF encoding VP2 is near the 3'-end of 
the virus genome and is conserved among caliciviruses [5] . 

The genus Vesivirus currendy contains two approved 
species. Vesicular exanthema of swine virus (VESV) and 
Feline calicivirus (FCV), and a diverse group of unassigned, 
phylogenetically-related viruses [6]. There are several well- 
recognized animal pathogens among vesiviruses that have 
been associated with a variety of disease conditions. These 
include diarrheal disease in dogs [7], respiratory illness, ves- 
icular lesions, and epidemic hemorrhagic fever in cats [8,9], 
and vesicular lesions in several other host species including 
swine, pinnipeds and humans [10-12]. The prototype virus 
of the genus Vesivirus, VESV, was originally isolated from 
pigs with clinical signs compatible with those caused by 
infection with foot-and-mouth disease virus (FMDV) 



o 



© 2012 Sandoval-Jaime et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
BIOIVIGCI Central creative commons Attribution License {http://creativecommons.0rg/licenses/by/2.O), v^/hich permits unrestricted use, 
distribution, and reproduction in any medium, provided the original w/ork is properly cited. 



Sandoval-Jaime et al. Virology Journal 2012, 9:297 
http://www.virologyj.eom/content/9/1/297 



Page 2 of 1 1 



[11]. In 1972, a virus with morphological and biochem- 
ical characteristics indistinguishable from those of VESV 
was isolated on San Miguel Island, California from sea 
lions and named San Miguel Sea lion virus (SMSV) 
[13,14]. When the experimental infection of pigs with 
SMSV resulted in a vesicular disease clinically mimicking 
FMDV and VESV infection, marine animals including 
ocean fish were retrospectively implicated to be the 
original source of VESV outbreaks [13,14]. Since then, 
VESV SMSV and other related caliciviruses have fre- 
quently been designated as the "marine vesiviruses." In 
contrast to FCV, which are considered to have 
restricted host specificity to cats of the family Felidae, 
the marine vesiviruses have been described as having 
an unusually broad host range [15-18]. 

In 1978-79, sixteen new vesivirus strains were isolated 
from four poikilothermic species (Aruba Island rattlesnake, 
Crotalus unicolor, Rock rattlesnake, Crotalus lepidus, 
Eyelash viper Bothrops schlegeli, Bell's horned frog Cere- 
tophyrs ornate) in a California zoological collection. The 
sixteen viruses were antigenically related and were not 
neutralized by the available VESV-like reference sera 
[19]. The new viruses were proposed as members of a 
new reptilian caliciviruses (RCV) Crotalus 1 (Crol) sero- 
type [19]. Sequence analysis of a 453 nt region of the Crol 
polymerase gene provided additional evidence for a new 
vesivirus group [20]. The Crol serotype did not appear to 
be restricted geographically or temporally, or limited to 
reptile and amphibian hosts. In 1986-7, vesiviruses neutra- 
lized by the Crol typing serum were isolated from samples 
collected from three different marine mammals species 
(Eumetopias jubatus, Zalophus californianus californianus, 
and Callorhinus ursinus) along the coast of Oregon and 
California states [16]. 



In this study, we determined the full-length genome 
sequence of Crol strain 7800321, isolated from the intes- 
tine of a Rock rattlesnake {Crotalus lepidus) housed in 
the San Diego Zoo in 1978. Comparison of the genome 
sequence of 7800321 with those available in GenBank 
database shows that this Crol virus represents a genetic- 
ally distinct vesivirus strain within the species VESV. 



Results and discussion 

Crol virus stocks obtained from snake samples collected 
during the original outbreak in the San Diego Zoo were 
examined for their ability to grow in cell culture following 
decades of storage at -80°C. In a preliminary screening, 
virus sample #7800321 obtained from a Rock rattlesnake 
{Crotalus lepidus) showed efficient growth in Vero cells 
(CCL-81, ATCC, Manassas, VA). Strong cytopathic effect 
(CPE) was observed in the virus infected cell monolayer at 
24-48 hours post infection (hpi). Plaque titration with an 
agarose overlay resulted in the formation of detectable pla- 
ques at 48 hpi (Figure lA). The efficiency of virus replica- 
tion was examined in a multiple-cycle growth curve time 
course analysis. Inoculation of Vero cells at multiplicity of 
infection (MOI)=0.01 consistently produced titers of -10*- 
10^ pfu/ml by 24 hours (Figure IB). To verify virus identity, 
four individual virus plaques were collected in agarose plugs 
and virus RNA was extracted from each of them using the 
RNeasy Kit (Qiagen, Valencia, CA). Purified RNA was 
employed for the RT-PCR amplification of the virus subge- 
nomic region followed by direct sequencing of the ORF2. 
Comparison of the ORF2 nucleotide sequences with those 
in GenBank (AY772540 and AY772541) confirmed the 
identity of the 7800321 strain as Crol. No evidence of gen- 
etic variation was observed among plaque-purified viruses. 



B. 




1010 




0 8 16 24 

hours 

Figure 1 In vitro growth characterization of the 780032! strain isolated from Crotalus lepidus. A) The growth of the 7800321 strain results 
in lysis of infected Vero cells and in the subsequent formation of plaques in a cell monolayer overlayed with 1% agarose containing growth 
medium. Plaques can be visualized with crystal violet staining after cells are fixed with formaldehyde. B) The virus titers of the 7800321 strain 
produced by the Vero cells. The cell monolayers were infected with MOI=0.01 and virus titer was measured at different time points by plaque 
assay. The titers shown are expressed as the mean from two independent experiments performed in duplicate. 
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Re-amplified 7800321 stock was employed for virus 
RNA isolation, cDNA amplification and direct sequen- 
cing analysis. To determine the full-length genome se- 
quence, several overlapping cDNA fragments of the 
virus genome were RT-PCR amplified, and sequenced 
on both strands using a primer-walking technique. The 
complete genome sequence of the 7800321 strain con- 
tained 8296 nt, excluding the poly (A) tail. Like other 
vesiviruses, it was comprised of three predicted open 
reading frames (ORFl, 2, 3) bordered by short 5'- and 
3'-end nontranslated (NTR) sequences (Figure 2). The 
5'-end NTR sequence of the 7800321-strain genome was 
19 nt in length and was highly conserved among VESV- 
related vesiviruses (VESV A48, PANl, RaV, WCV, 
SMSVl, v810 and vl415). The 3'-end NTR was 164 nt 
in length, and showed a higher degree of nucleotide vari- 
ation when compared to that of other vesiviruses. The 
lengths of the ORFl, ORF2 and ORF3 were found to 
be 5643 nt, 2121 nt and 348 nt, respectively. The stop 
codon of ORFl (in bold) and the start codon of ORF2 
(underlined) were separated by five nucleotides 
(UAGCCAUU AUG ), while ORF2 and ORF3 overlapped 
by four nucleotides ( AUG A) similar to other vesi- 
viruses. The 3'-end of the 7800321 ORFl (nt 5648-5660) 
contained a conserved sequence motif that showed a 
high level of identity (10 out of 13 nt) with the first 13 nt 
of the virus genome. In addition, this motif was highly 
similar (9 out of 13 nt) to the 5'-end of the FCV subge- 
nomic RNA that was mapped by primer extension 
analysis [21,22]. Taken together, these data supported 
the predicted start of the 7800321 subgenomic RNA as 
position 5648 of the genome. 



The vesivirus ORFl encodes a nonstructural polyprotein 
that undergoes co-translational proteolytic processing 
during virus replication [23]. A cascade of proteolytic 
events mediated by the virus-encoded proteinase gives rise 
to the virus mature nonstructural proteins and their inter- 
mediate forms. The Crol ORFl encodes an 1880 amino 
acid nonstructural polyprotein with a predicted molecular 
mass of 209 kDa. Motif scanning [24,25] and comparative 
sequence analysis of the Crol polyprotein revealed the 
presence of the characteristic NTPase ( GppgcGKT ), 
3C-like proteinase (^^"^GDCG""'') and 3D-like RNA- 
dependent RNA polymerase (^^^^GLPSG^^^^ and 
i70iY(t;]-)£)i704^ catalytic motifs in an order conserved 
among caliciviruses. Alignment of the Crol ORFl poly- 
protein sequence with that of another vesivirus, FCV, that 
had an experimentally established cleavage map [23] 
allowed prediction of the putative cleavage sites and, 
correspondingly, the sizes of mature nonstructural proteins. 
The predicted Crol polyprotein cleavage sites were con- 
served among VESV-related vesiviruses and, with the ex- 
ception of cleavage site between NS3 and NS4 proteins, all 
carried a glutamic residue in the PI position and alanine, 
serine or glycine residues at the position PI' (Figure 2). Of 
interest, the NS3-NS4 scissile bond was predicted between 
glutamine and alanine residues. The order and sizes of the 
proteins defined by these cleavage sites were as follows: 
16.2 kDa NSl - 32.2 kDa NS2 - 39.6 kDa NSS^™""' - 
31.2 kDa NS4 - 13.5 kDa NS5^s _ 76,5 NS6-7''"''°'. 

The vesivirus ORF2 encodes a precursor of the virus 
capsid protein that is processed by cleavage during cap- 
sid protein maturation. Processing of the FCV capsid 
precursor at an E'^^'^-A^^^ dipeptide by virus NS6-7'''^°''°' 
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Figure 2 Schematic representation of the 7800321 strain genome organization. The 7800321 genome is comprised of three major ORFs. 
ORFl encodes a nonstructural polyprotein, 0RF2 - a precursor of the virus capsid protein, VPl, and ORFS - a minor capsid protein, VP2. Mature 
virus proteins are shown as rectangular boxes with calculated protein molecular masses indicated below. The putative borders of the virus 
mature nonstructural proteins are indicated with predicted cleavage sites and arrows. The E^^VS site that is likely cleaved by the virus-encoded 
protease during capsid protein maturation is shown for the ORF2-encoded precursor. The 7800321 predicted cleavage sites in the ORFl were 
based on an alignment of VESV-like ORFl polyprotein sequences with the established map of FCV [23]. 
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results in removal of an 124 aa-long N-terminal leader 
sequence (LC) and release of the 59-60 kDa mature 
form of the virus capsid protein, VPl [26]. The Crol 
ORF2 sequence is predicted to encode a precursor 
protein with a molecular mass of -78 kDa. Alignment 
with the FCV capsid precursor sequence showed the 
presence of a putative scissile bond between ^^^E and 
'^^^S (Figure 2). It is likely that, similar to FCV, this 
cleavage site defines the border between the virus capsid 
protein leader sequence and mature VPl. Of interest, the 
site is conserved in the corresponding sequences of the 
recently characterized v810 and vl415 vesivirus strains 
isolated from Steller sea lions [27]. The expression of the 
vSlO and vl415 ORF2 sequences that were N-terminally 
truncated at this site resulted in production and self- 
assembly of virus-like particles morphologically and anti- 
genically similar to virions [28] . 

The small ORF3 located at the 3'-end of the vesivirus 
genome encodes a minor component of virus capsids, the 
VP2 protein [29]. The protein is expressed from the virus 
bicistronic subgenomic RNA and has been shown to be es- 
sential for the formation of infectious vesivirus virions [30]. 
The substitution of a C for a T at nucleotide 8115 con- 
verted a stop codon (TAG) to a glutamine codon (CAG) 
resulting in a five codon-extension of the ORF3 compared 
to other VESV-related vesiviruses (data not shown). The 3'- 
end of the calicivirus RNA genome has been shown to play 
a crucial role in the initiation of virus replication [31,32]. 
Nevertheless, modifications introduced into the genomes 
of feline calicivirus and murine norovirus showed that this 
region could tolerate a number of sequence changes while 
supporting virus replication [30,32]. 

To elucidate the genetic and phylogenetic relationships 
of the Crol virus with other vesiviruses, multiple se- 
quence alignments of the representative vesivirus genome 
sequences were generated using the ClustalW algorithm 
in the Mega5 software package [33]. Comparative genetic 
distance analysis of the vesivirus sequences revealed the 
presence of three genetically distinct virus groups, provi- 
sionally called here as VESV-like (or marine vesiviruses), 
FCV-like and CaCV-like viruses. The newly characterized 
Crol virus clustered with the VESV-related group, where 
differences in sequence identity did not exceed 23%. Pair- 
wise comparison of the Crol genome sequence with 
those of VESV-Uke viruses showed that the genome of 
the most distant VESV-like virus, walrus calicivirus, dif- 
fered from that of Crol by 21%. The genetic distances 
observed between genomes of Crol and FCV-like, and 
CaCV-lUce viruses were significantly higher at 48% and 
50%, respectively (Table 1). Of interest, FCV and CaCV 
were most distant within the genus, with differences in 
genome sequences reaching 56% (data not shown). 

The genome variability of viruses within the VESV-like 
group was analyzed using the Plotcon algorithm in the 



EMBOSS software package [34]. A similarity profile 
generated for the multiple sequence alignment of 
VESV-like viruses revealed the presence of conserved 
and variable regions across the virus genome (Figure 3). 
The relatively well conserved sequences included virus 
nonstructural NSl and NS4 genes, and the region cor- 
responding to the ORF1-ORF2 junction known to contain 
a putative transcription start of the virus subgenomic 
RNA. The VESV-like NSl and NS4 genes shared an 
increased level of average nucleotide sequence identity 
(88.6%) compared to that of the entire ORFl (84.6%). As 
expected, the most variable part of the virus genome was 
observed in the ORF2 region that encodes the P domain 
of the virus capsid protein. The marked variation (35.4% 
nucleotide difference) was consistent with structural and 
functional studies that have shown that this region of the 
VPl is exposed on the surface of the virion and plays a 
role in receptor binding and immune selection in the host 
[35]. Of interest, similarity plots for the deduced amino 
acid sequences encoded by ORFl, 2 and 3, had profiles 
similar to those of the corresponding nucleotide regions 
of the virus genome (data not shown). 

The profiles of the genome plots generated in a compari- 
son of VESV-, FCV- and CaCV-like vesivirus sequences as 
well as that for VESV-like sequences alone were similar. 
Consistent with the generally higher level of genetic diver- 
gence existing between vesiviruses from different groups, 
the similarity indexes calculated for the all-vesivirus se- 
quence alignments were significandy lower. Interestingly, 
the comparative analysis showed that the NSl and part of 
the NS4 sequence were not conserved when compared 
to the corresponding sequences of the FCV-like and 
CaCV-like viruses (Figure 3). Furthermore, predicted vesi- 
virus NSl gene sequences varied in length from 138 to 
534 nt between these groups of viruses. For the Crol, 
the nucleotide identity of the NSl gene (predicted to be 
148 aa long) was 29.7-38.5% with the corresponding 
sequences of CaCV-like viruses, and 51.6-52.2% with 
those of FCV-like viruses. The level of similarity for the 
NSl deduced amino acid sequences was calculated to 
be 19.6-22.2% and 8.8% between Crol and CaCV-like 
and between Crol and FCV-like viruses, respectively. 
The function of the calicivirus NSl protein is un- 
known and cannot be predicted based on the protein 
sequence since it does not show significant homology 
with any established functional sequence motifs. More- 
over, the level of similarity of NSl amino acid 
sequences between different calicivirus genera does 
not exceed that of random sequences, precluding identi- 
fication of conserved sequence motifs. 

The observed level of the NS4 sequence identity was 
significantly lower when the Crol sequence was aligned 
to those of FCV- and CaCV-like viruses. For the latter 
two, it ranged from 55.2 (Crol vs FCV) to 56.4% (Crol 



Table 1 Percent nucleotide and amino acid identity of Cro1 strain 7800321 with other versiviruses 
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Relatedness of amino acid sequences are shown as percent identity/percent similarity. 
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Figure 3 Schematic representation of the vesivirus genome organization and the sequence conservation along the alignment of the 
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vs CaCV). In contrast, the lowest level of sequence iden- 
tity among VESV-like virus NS4 genes was 86.4%. The 
similarity profile showed that the most variable 
sequences of the NS4 gene were located near the 5'-end 
(Figure 3). Similar to NSl, the NS4 gene is not con- 
served among viruses from different Caliciviridae genera 
and represents the second most variable region in the 
nonstructural ORF of their genomes. Nevertheless, all 
calicivirus NS4 proteins share a conserved structural 
feature, which is the presence of a hydrophobic domain. 
Of interest, a cluster of hydrophobic amino acid residues 
is located near the C-terminus of the Crol NS4 (Mac- 
Vector's Protein Analysis Toolbox), with amino acids 
249-271 predicted to form a membrane-associated helix 
(TMpred server). Consistent with the putative role of this 
domain in membrane interactions, biochemical studies 
showed that transiently expressed NS4 behaved as an 
integral membrane protein [36]. In addition, different 
forms of the NS4 were found to localize to the 
membrane-associated virus replication complexes in 
calicivirus infected cells [36,37]. The significant sequence 
diversity suggests that NS4 might play an important role 
in determining the specificity of protein-membrane inter- 
actions in the host cell. For example, the subcellular 



localization of the norovirus NS4 is determined by an ER 
export signal motif (MERES) conserved only among the 
noroviruses [38]. Moreover, the presence of the MERES 
motif was shown to be critical to the NS4 antagonist role 
in ER/Golgi trafficking [38,39]. Of interest, computa- 
tional analysis showed that the MERES motif was not 
present in the Crol sequence. In addition, scanning of 
the Crol NS4 sequence with software designed to 
identify putative signal and subcellular localization motifs 
(see Materials and Methods) found no known targeting 
sequences. The presence of such signals in the Crol NS4 
protein remains to be established. 

Another region of marked sequence variation was 
observed downstream from the ORF1-ORF2 junction. 
Plotcon analysis of the vesivirus sequences revealed a 
low level of nucleotide identity among the virus LC 
genes, with the lowest identity (30.8%) between the Crol 
and FCV65 viruses. Accordingly, a lower level of similarity 
was found for the compared deduced amino acid 
sequences of this protein, 30.1-30.7% for FCV-Crol and 
39.4-43.1% for CaCV-Crol pairs. The function of this pro- 
tein remains unknown; however, cleavage of the LC from 
the capsid precursor molecule was crucial for production 
of infectious virus particles [26] . 
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To investigate the phylogenetic relationship of the 
Crol virus with other vesiviruses, a phylogenetic tree 
was inferred from multiple alignments of representative 
sequences using the Bayesian method. Figure 4A shows 
a consensus tree produced by MrBayes3.1.2 for the set 
of full-length genome sequences. A similar approach 
was employed to generate additional phylogenetic trees 
for the vesivirus ORFl and subgenomic RNA sequences 
(Figure 4B and 4C). Analysis of the vesivirus genomic 
tree revealed the presence of three phylogenetic groups. 
The first group included FCV strains; the second group 
consisted of viruses related to the CaCV strain, and the 
third combined together viruses closely related to the 
genus prototype strain, VESV A48. A similar topology 
was observed for the phylogenetic trees inferred for the 
ORFl and subgenomic RNA sequences. In all trees, the 
presence of three major clusters was strongly supported 
by a high posterior probability value for each clade. 

Consistent with genetic distance analysis, the genomic 
sequence of the Crol 7800321 strain clustered together 
with those of the VESV-like viruses, making it a member 
of the marine vesivirus group (Figure 4). Within this 
cluster, the inferred phylogenetic relationships between 
RaV, SMSVl, VESV A48, PANl, WCV and Steller sea 
lion v810 and vl415 strains were found to be similar to 
those reported earlier [17,27,40]. Of interest, when the 
phylogenetic tree was inferred for the virus subgenomic 
RNA sequences, the 7800321 strain clustered together 
with the v810 strain. Similar clustering was observed 



when a set of the capsid sequences was extended to in- 
clude two additional sequences available in GenBank, 
SMSV4 and SMSVl 7 (accession numbers M87482 and 
U52005, respectively). In the extended capsid tree, 
SMSV4 formed a separate group with VESV A48, while 
SMSV17 clustered with the 7800321, v810 and vl415 
strains. However, the latter cluster showed a lack of 
strong statistical support, with the posterior probability 
value reaching only 0.68 (data not shown). Of note, the 
ORFl and subgenomic sequences of the SMSVl strain 
showed inconsistent grouping with those of RaV and 
WCV. The observed incongruence of the phylogenetic 
clustering of these viruses suggested a possible recom- 
binant origin of the SMSVl strain. Similarity plot analysis 
of the SMSVl and other vesivirus genome sequences 
demonstrated an increased level of nucleotide sequence 
identity between the subgenomic regions of SMSVl and 
RaV/WCV strains with a predicted putative recombin- 
ation site near the junction of the SMSVl ORFl and 
ORF2 sequences (data not shown). However, our prelim- 
inary analysis could not identify the second parental 
strain or its probable lineage (data not shown). Confirm- 
ation of a possible recombination event will require 
further investigation. 

More than 40 serotypes of vesiviruses have been 
identified using serum-neutralization tests [41]. The 
virus capsid protein VPl has been shown to encode 
major antigenic epitopes recognized by antibodies in 
polyclonal sera of infected animals. Correspondingly, 
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the sequence variability of this protein is thought to 
provide the molecular basis for serotypic diversity of 
the vesiviruses. The highest similarity of amino acid 
sequence (91.4%) for the Crol VPl protein was 
observed with the VPl protein of a recently described 
v810 strain [28]. Interestingly, most (72 out of 79) of 
the amino acid changes in the VPl sequence of these 
strains mapped to the C-terminal part of the protein, a 
region of the vesivirus genome with the highest level of 
variability (Figure 3). Evolutionary conservation scores 
generated from the corresponding multiple sequence 
alignment were mapped onto the structure of the 
SMSV4 VPl protein with the aid of the ConSurf server 
software [42]. The X-ray structure of the SMSV4 VPl 
(resolved to 3.2-A) showed that the vesivirus capsid 
protein shared a domain organization similar to that of 
human noroviruses. The protein contained a short N- 
terminal arm (aa 10-48) and two domains, shell (S) 



and protruding (P), formed by its N-terminal (aa 49-209) 
and C-terminal (aa 210-554) parts, respectively. Based on 
the structure, the P domain could be further divided into 
two subdomains, PI and P2, with the latter formed by 
the most distally located part of the protein [35]. The 
ConSurf analysis of the sequence variability coupled 
with visualization of the data in Chimera [43] showed 
that the most divergent parts of the P domain were loca- 
lized on the surface of the virus capsid and were repre- 
sented mainly by the sequence of the P2 subdomain 
(Figure 5 A and 5B). It was noteworthy that 54 out of 72 
amino acid differences observed between the VPl P 
domains of the v810 and 7800321 strains were mapped 
to the P2 region (Figure 5A). Structurally, the P2 sites 
that accumulated the majority of the sequence differ- 
ences between these two viruses were located in five 
loops exposed on the surface of the virus capsid 
(Figure 5 A and 5B). Surface localization of these loops 
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and their flexibility are thought to play an important 
role in defining antigenic characteristics of the virus 
capsid and in its interactions with host cell receptors. 

Our finding that Crol strain 7800321 (isolated from 
Rock rattlesnake) shares strong sequence identity with 
the VESV-like marine vesiviruses is consistent with 
reports that the Crol serotype has become established 
in both aquatic and terrestrial hosts [16,19]. The original 
tissue samples fi"om reptiles in the San Diego Zoo Crol 
outbreak of 1978 are no longer available for analysis, so it 
is impossible to investigate the outbreak retrospectively. 
Although the source of this Crol virus remains un- 
known, the continued genetic characterization of known 
viral isolates in concert with the increasing use of deep 
sequencing techniques for virus discovery may help track 
the origin and spread of this and other caliciviruses in 
nature. 

Materials and methods 

Sample isolation and virus amplification 

Reptile tissue samples (small intestine, liver, kidney) were 
collected in 1978 from hatchling and breeding snakes in 
the San Diego Zoo. Snakes were experiencing high mor- 
tality rates thought to be associated with enteritis and 
hepatitis of an unknown cause [19]. The tissue samples 
were homogenized and clarified by low-speed centrifuga- 
tion. To amplify the virus, the corresponding superna- 
tants were filtered through 0.22 i^m polysulfone filters, 
and the resulting filtrates were added to Vero cell mono- 
layers maintained in Eagle's minimal essential medium 
that was supplemented with penicillin (200 units/ml), 
streptomycin (100 j^g/ml), L-glutamine (2 mM) and 10% 
heat-inactivated fetal bovine serum. The inoculated cells 
were monitored until the appearance of visible CPE. 
When CPE in monolayers exceeded 90%, virus stocks 
were generated by collection of growth media and low 
speed clarification of supernatants. After viruses were 
plaque-purified three times, the amplification procedure 
was repeated and the clarified supernatants were ali- 
quoted and stored at -80°C. 

RNA extraction and RT-PCR sequencing 

Viral RNA was extracted from virus stock samples with 
the RNeasy Mini Kit (Qiagen). The extracted RNA was 
used as a template for reverse transcription (RT) and 
PCR amplification of cDNA fragments. The RT-PCR 
reactions were performed using One-Step RT-PCR Kit 
(Invitrogen, Carlsbad, CA) and vesivirus genome- 
specific primer pairs. Briefly, following the initial RT 
reaction, 30 min at 45°C, a denaturation step at 94°C 
was performed for 2 min, followed by 40 cycles of 15 s 
at 94°C, 30 s at 50°C, and 3 min at 68''C. Amphfied 
cDNA fragments were resolved by electrophoresis in 
1.0% agarose gels containing ethidium bromide. The 



corresponding DNA bands were visualized with UV 
light and excised from the gel. The DNA was extracted 
with QIAquick Gel Extraction Kit (Qiagen) and sub- 
jected to nucleotide sequencing using the Big Dye Ter- 
minator v3.1 Cycle Sequencing Ready Reaction Kit and 
an automated sequencer, ABI 3100 (Applied Biosys- 
tems, Carlsbad, CA). Sequences of the primers that 
were employed for cDNA fragment amplification are 
given in Additional file 1: Table SI and sequences of 
the primers used for genomic sequencing are available 
upon request. The sequences of the 5'-end regions of 
the virus genome were determined using 573 RACE 
Kit 2"'' Generation (Roche). The corresponding cDNA 
fragments were synthesized and amplified using 5'-end 
sequence-specific and anchor (5'-GACCACGCGT- 
ATCGATGTCGACTTTTTTTTTTTTTTTTV-3') primers 
according to protocol provided by the manufacturer. 
Following agarose-gel purification, fragments were sub- 
jected to direct sequence analysis as described above. 
Virus genome sequences were assembled using the 
Sequencher 4.9 program (Genecodes, Ann Arbor, MI). 

Computer sequence analysis 

Derived Crol nucleotide and amino acid sequences were 
analyzed, aligned and compared with vesivirus sequences 
available from the GenBank database using MacVector 
(MacVector, Inc., Gary, NC), EMBOSS [34] and Mega5 
[33] software packages. The GenBank accession numbers 
for complete vesivirus genome sequences were: Canine 
calicivirus (CaCV), AB070225; Vesicular exanthema of 
swine virus (VESV A48), U76874; San Miguel Sea Lion 
virus-1 (SMSVl), AF181081; Primate calicivirus 1 (Panl), 
AF091736; Walrus calicivirus (WCV), AF321298; Rabbit 
vesivirus (RaV), AJ866991; Feline calicivirus F9 strain 
(FCVF9), M86379; Feline calicivirus F4 strain (FCVF4), 
D31836; Feline calicivirus F65 strain (FCVF65), AF109465; 
Feline calicivirus Urbana strain (FCVURB), L40021; Steller 
sea lion vesivirus-v810 (v810), EF193004; Steller sea lion 
vesivirus-vl415 (vl415), EF195384; Calicivirus isolate AU- 
ston 2008/US (ACV8), GQ475302; Calicivirus isolate Geel 
2008/Belgium (GCV8), GQ475303; Calicivirus isolate AU- 
ston 2009 /US (ACV9), GQ475301). The genome sequence 
of the 7800321 strain was submitted to GenBank and 
assigned accession number JX047864. 

A plot of average similarity for each set of the aligned 
vesivirus sequences was generated using the Plotcon 
program from the EMBOSS software package (http:// 
emboss.open-bio.org/ rel/ rel6/ apps/plotcon.html). 

Scanning for the protein domains and motifs in the 
predicted protein sequences were performed using Inter- 
ProScan (http://www.ebi.ac.uk/Tools/pfa/iprscan). Signal 
peptide and subcellular localization motifs predictions 
were performed using iPSORT (http://hc.ims.u-tokyo. 
ac.jp/iPSORT), PSORTII (http://psort.nibb.ac.jp/form2. 
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html), SOSUIsignal (http://bp.nuap.nagoya-u.ac.jp/sosui/ 
sosuisignal/sosuisignal_submit.html), SIG-Pred (http:// 
bmbpcu36.leeds.ac.uk/prot_analysis/Signal.html), Golgi 
Predictor (http://ccb.imb.uq.edu.au/golgi/golgi_predictor. 
shtml), PTSl predictor (http://mendel.imp.ac.at/mendelisp/ 
sat/ptsl/PTSlpredictor.jsp), Predotar (http://urgi.versailles. 
inra.fr/predotar/predotar.html), and SignallP (http:// 
www.cbs.dtu.dk/services/SignalP). Prediction of the mem- 
brane associated domains was performed using 
TMpred Server (www.ch.embnet.org/software/TMPRED_ 
form.html). 

Bayesian inference of phylogeny was carried out with 
Mr.Bayes 3.1.2 software [44]. The analysis was per- 
formed using Markov chain Monte Carlo sampling 
under a general time-reversible model of nucleotide 
substitution with a gamma distribution of rates and a 
proportion of invariant sites. The search was run for 
one million generations with tree sampling occurring 
every 100th generation. Bayesian posterior probabilities 
and tree topologies were calculated from the consensus 
of collected tree samples after excluding the first 25% 
trees as burn-in. The analyses were performed two 
times for every set of the sequences. The resulting trees 
were visualized with FigTreel.3.1 (http://tree.bio.ed.ac. 
uk/software/figtree/). 

Generation of multiple-cycle growth curve 

Vero cell monolayers in sbc-well plates were infected with 
the 7800321 strain at MOI=0.01, and virus titers were 
measured at different time points post infection. Briefly, 
virus stocks were diluted in a supplemented medium (see 
above), added to cells and allowed to adsorb for one hour 
at 37 °C. After the incubation, the inocula were removed, 
the infected cells were washed and fresh growth medium 
was added to each well. At 0, 8, 16, and 24 hpi, the plates 
were frozen and stored at -80 °C. All collected samples 
were freeze-thawed three times, and cell debris was 
removed by low speed centrifugation. Aliquots of the 
clarified supernatants were used to determine virus titers 
at each time point by plaque assay [45]. 

Additional file 



Additional file 1: Table SI. PCR primers used in amplification of the 
Crol genome cDNA fragments. 
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